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I. INTRODUCTION 

In this paper, we describe classical field theories by means of a new variational principle, 
referred to as the Hamilton-Pontryagin variational principle and we show that the resulting 
field equations can be described in terms of multi-Dirac structures, a field theoretic extension 
of the concept of Dirac structures introduced by Courant in Ref. 14. To set the stage for 
the paper, wc begin with reviewing the Hamilton-Pontryagin principle and Dirac structures 
in the context of mechanical systems. 



A. The Hamilton-Pontryagin Variational Principle in Mechanics 

Let be a Lagrangian. The Hamilton-Pontryagin principle is a variational princi- 

ple in which the position coordinates and the velocity coordinates are treated indepen- 
dently and the relation = f * is imposed as a constraint by means of a Lagrange multiplier 
Pi- This leads to an action functional of the form 

S{q/o,p) = f ^ {L{q\t),v\t)) + {p,{t) , q\t) - v\t)) ) dt 

(1) 

{{Piit),m) - E{q'{t),v'it),p,it)))dt, 



to 

where E{q^, v^,Pi) :— {pi, v^) —L{q^, v^) is the generalized energy associated to the Lagrangian. 
By taking arbitrary variations with respect to q\ and Pi, we obtain the Euler-Lagrange 
equations in implicit form: 

Pi ^Q-, P^=Q^v q ^v. (2) 

Variational principles in which the position and velocity coordinates are varied indepen- 
dently have a long history: the action functional (1) in particular first appears in the work 
of Livens in Ref. 29 and was rediscovered in Refs. 20 and 35, among others. The name 
"Hamilton-Pontryagin principle" was coined by Yoshimura and Marsden in Ref. 42 because 
of the similarity with the Pontryagin maximum principle in optimal control theory. 

It was shown in Ref. 42 that the Hamilton-Pontryagin principle and the associated implicit 
equations of motion can be described in an intrinsic way by means of Dirac structures. A 
discrete version of the Hamilton-Pontryagin principle was developed in Refs. 4 and 27 with 
a view towards the design of accurate variational integrators. 
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B. Dirac Structures and Lagrange-Dirac Systems in Mechanics 

The notion of Dirac structures was originally developed by Courant and Weinstein in 
Refs. 13 and 14 and Dorfman in Ref. 15 as a unification of (pre-)symplectic and Poisson 
structures. It was soon realized that Dirac structures play an important role in mechanics. 
In particular, it was shown that interconnected systems, such as LC circuits, and nonholo- 
nomic systems can be effectively formulated in the context of implicit Hamiltonian systems 
(Refs. 13-15, and 33). On the Lagrangian side, it was shown in Ref. 41 that Dirac structures 
induced from distributions on configuration manifolds naturally yield a notion of implicit 
Lagrangian systems, allowing for the description of mechanical systems with degenerate 
Lagrangians and nontrivial constraint distributions. 

In order to show the link with the equations (2) obtained from the Hamilton-Pontryagin 
principle, let Aq be a distribution on a manifold Q. Consider the canonical symplectic form 
Q on T*Q and define Qm as the puUback of fl to the Pontryagin bundle M :— TQ © T*Q. 
The set D, given by 

D = {{X, a)eTM® T*M \ i^^^M - a e A^, X e Am}, 

where Am := ^^qmIAq) (with tvq^m : M — > Q is the Pontryagin bundle projection), is a 
Dirac structure (see Refs. 41 and 42); namely, D is maximally isotropic with respect to the 
following pairing on TM © T*M: 

{{{X,a),{Y,/3)))_^l{ix/3 + iya). (3) 

The implicit Euler- Lagrange equations are given by {X, dE) e D. It can be shown 
that these equations are given in coordinates by 

(9/ (9/ 

Pi=9^V Pr=Q^i+^aAt{q), (i = ^, and ^^^K = 0; (4) 

see Ref. 12. Here, the forms A^{(i)d(t form a basis for the annihilator A°, and the are 
Lagrange multiphers. The equations (4) are readily seen to be the nonholonomic Euler- 
Lagrange equations in implicit form, and when the constraints vanish, = 0, we recover 
the imphcit Euler-Lagrange equations (2). 
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C. The Hamilton-Pontryagin Principle for Field Theories 

For classical field theories, the Lagrangian L{x^, y^, v^) depends not only on the space- 
time coordinates and the fields y^, but also on the spacetime velocities v^, which repre- 
sent the derivatives of the fields y"^ with respect to x^^. By comparison with (1), a candidate 
Hamilton-Pontryagin principle may then be given in coordinates by 



r+'x. (5) 



By taking variations with respect to y^, and p'\ (where 5y^ is required to vanish on the 
boundary dU of the integration domain), we obtain the following implicit equations: 



A 



which clearly generalizes (2) to the case of field theories. In section III, we provide an 
intrinsic version of this variational principle, using the geometry of the first jet bundle and 
its dual as our starting point. In contrast to mechanics, the jet bundle is an affine bundle 
and this makes the definition of the duality pairing used implicitly in (5) somewhat more 
complicated. 



D. Multi-Dirac Structures 



In section IV, we describe the implicit field equations (6) in terms of multi-Dirac struc- 
tures. These geometric structures were introduced in Ref. 38 as a graded version of the 
standard concept of Dirac structures, and their relevance for field theory lies in the fact that 
they generalize the concept of the "graph of a multisymplectic structure." 

To make this more precise, we recall that a multisymplectic structure is a form Q of 
degree at least two, which is closed and non-degenerate. When the degree of Q is two, Q 
is a symplectic form, while if ^2 is a top form, then it is necessarily a volume form. The 
literature on multisymplectic field theories is by now very extensive, but for fundamental 
aspects we refer to Refs. 2, 5, 8, 21, 22, and 24 and the references therein. 

Just as the graph of a symplectic form is a Dirac structure, the graph of a multisymplectic 
form Q on a manifold M turns out to be an example of what we have called a multi-Dirac 
structure in Ref. 38. Here, the graph of Q is defined as follows. Let A; -|- 1 be the degree of Q 
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and consider the mapping whereby an Z-multivector field Xi (with 1 < Z < A;) is contracted 
with Q: 

The graph of this mapping determines a submanifold A of /\\TM) x m /\''^^~\T* M) , and 
by considering the direct sum 

D - © • • • © Dfc (7) 

of all such subbundles, we obtain a multi-Dirac structure. This structure is maximally 
isotropic under a graded antisymmetric version of the pairing (3), and multi-Dirac structures 
can be considered more generally as graded versions of standard Dirac structures. 

The link between multi-Dirac structures and the implicit field equations is provided 
by theorem III. 4, wherein we show that the field equations obtained from the Hamilton- 
Pontryagin principle can be written as Ix^m = (— l)"^^d£'. In multi-Dirac form this be- 
comes 

{X, i-ir+^dE) e (8) 

where Dn+i is the component of the highest degree in (7). More generally, we arrive in 
definition IV. 5 at the concept of a Lagrange- Dirac field theory, which is a triple {X,E, 
satisfying (8), where D^+i belongs to a multi-Dirac structure D which does not necessarily 
come from a multisymplectic form. 

E. Outline of the Paper 

After introducing the Hamilton-Pontryagin principle for field theories in section HI, we 
discuss the link with multi-Dirac structures in section IV. We finish the paper in section V 
with a number of examples of implicit field theories and multi-Dirac structures, taken from 
mechanics, electromagnetism and elasticity. Among others, we show that the Hu-Washizu 
principle from linear elastostatics is a particular example of the Hamilton-Pontryagin prin- 
ciple. 

The focus throughout this paper is on the Hamilton-Pontryagin principle and the ap- 
plication of multi-Dirac structures to field theories. In the companion paper Ref. 38, we 
develop the algebraic and geometric properties of multi-Dirac structures. 

5 



The Hamilton-Pontryagin Principle and Multi-Dirac Structures 



II. THE GEOMETRY OF JET BUNDLES 

In this section, we provide a quick overview of the geometry of jet bundles for the treat- 
ment of Lagrangian field theories. Most of the material in this section is standard, and can 
be found in Refs. 22 and 32 and the references therein. 

A. Jet Bundles 

Let X be an oriented manifold with volume form r], which in many examples is space- 
time, and let ttxy Y ^ X he a. finite-dimensional fiber bundle which we call the covariant 
configuration bundle. The physical fields are sections of this bundle, which is the co- 
variant analogue of the configuration space in classical mechanics. For future reference, we 
assume that the dimension oi X is n + 1 and that ttxy is a bundle of rank N, so that 
dimy = n -I- -I- 1. Coordinates on X are denoted x'', /x = l,2,...,n-|-l, and fiber coordi- 
nates Y are denoted by y^, A = 1, ...,N so that a section : X ^ F of ttxy has coordinate 
representation = {x^,y^{x)). We will also assume that X is equipped with a fixed 
volume form r) given in adapted local coordinates by 77 = d^'^^x :— dx^ A • • • dx^'^^. 

The analogue in classical field theory of the tangent bundle in mechanics is the first 
jet bundle J^Y, which consists of equivalence classes of local sections of ttxf, where we 
say that two local sections 0i, 02 of Y are equivalent at x e X if their Taylor expansions 
around x agree to the first order. In other words, 0i and 02 are equivalent if 0i(x) = 02 (a^) 
and Tj;0i = T-E02. It follows that an equivalence class [0] of local sections can be identified 
with a linear map 7 : T^X TyY such that Tttxy 07 = Idy^x- As a consequence, J^Y 
is a fiber bundle over Y, where the projection TTy^jiy : J^Y — > y is defined as follows: 
let 7 : T^X — >■ TyY be an element of J^Y , then t^y,j^y{i) '■— V- Coordinates on J^Y are 
denoted by where the fiber coordinates represent the first-order derivatives 

of a section. They are defined by noting that any 7 G J^F is locally of the form 



The first jet bundle J^Y is an affine bundle over Y, with underlying vector bundle 
L{TX, VY) of linear maps from the tangent space TX into the vertical bundle VY , defined 
as 




(9) 



VyY ^{ve TyY I Tt:xy{v) = 0}, for yeY. 
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Given any section (p : X ^ Y oi ttxy, its tangent map T!r0 at a; e X is an element of 
JyY, where y — Thus, the map x i-> Tj;(f) defines a section of J^Y, where now J^Y 

is regarded as a bundle over X. This section is denoted j^c/) and is called the first jet 
prolongation of (p. In coordinates, j^(f) is given by 

where di, — d/dx^ . A section of the bundle J^Y — > X which is the first jet prolongation of 
a section : X — > y is said to be holonomic. 

B. Dual Jet Bundles 

Next, we consider the field-theoretic analogue of the cotangent bundle. We define the 
dual jet bundle J^Y* to be the vector bundle over Y whose fiber aX y E Y^ is the set of 
affine maps from J^Y to A^+^X, where A^~^^X denotes the bundle of {n + l)-forms on X. 
Note that since the space of affine maps from an affine space into a vector space forms a 
vector space, J^Y* is a vector bundle despite the fact that J^Y is only an affine bundle. A 
smooth section of J^Y* is an affine bundle map of J^Y to A"+^X. Any affine map from 
JyY to A^'+^X can locally be written as 

vi^{p + p'X)d-^'x, 

where := dx^ A dx'^ A • • • A dx""^^ is a coordinate representation of the volume form 77, 

so that coordinates on J^Y* are given by (x'*, 

Throughout this paper we will employ another useful description of J^Y*. Consider again 
the bundle A"+^y of (n + l)-forms on Y and let Z C A^+^y be the subbundle whose fiber 
over 7/ e F is given by 

Zy = {ze A'^+^Y I i^i^z = for all v,w E VyY}, 

where i„ denotes left interior multiplication by v. In other words, the elements of Z vanish 
when contracted with two or more vertical vectors. The bundle Z is canonically isomorphic 
to J^Y* as a vector bundle over Y; this can be easily understood from the fact that the 
elements of Z can locally be written as 

7 
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where d'^x^ :— J d^'^^x. 

Prom now on, we will silently identify J^Y* with Z. The duality pairing between elements 
7 e J,^ y and z e Zy^ can then be written as 

In coordinates, (z,^) — (p^^t"^ +p)d^'^^x. 

C. Canonical Multisymplectic Forms 

Analogous to the canonical symplectic forms on a cotangent bundle, there are canonical 
forms on J^Y*. Let us first define the canonical [n + l)-form 6 on J^Y* = Z hy 

0{z){ui, ...,Un+l) = z{T'Ky,j1y*{ui), TTTy^jiy* (lin+i)) ^^^^ 
= (7ry,Jiy*2;)(«l, ■■■,Un+l), 

where we have interpreted z & Z ^ J^Y* as before as an (n+l)-form on Y, and ui, Un+i £ 
TzZ. The canonical multisymplectic {n + 2) -form Q on J^Y* is now defined as 

n = -de. 

In coordinates, one has the following expression for ©: 

e = p^dy^ A d'^Xf, + p d^'+^x, 

while Q is locally given by 

n^dy^A dp'X A d"a;^ - dp A d"+^x. (11) 

It is easy to show (see Ref. 6) that Q is non-degenerate in the sense that, for all v e TZ, 
i^Q = implies that v — 0. Moreover, Q is trivially closed: dQ = 0. Porms which are both 
closed and non-degenerate are referred to as multisymplectic forms. If the degree of O is 
two, fl is just a symplectic form, while if Q is of maximal degree, ^2 is a volume form. More 
examples of multisymplectic manifolds can be found in Ref. 6. 

A form which is non-degenerate but not necessarily closed is referred to as an almost 
multisymplectic form. On the other hand, in the context of field theory, we will often 
refer to a closed but possibly degenerate form as a pre-multisymplectic form. 
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D. Lagrangian Densities and the Covariant Legendre Transformation 

A Lagrangian density is a smooth map C — Lrj: J^Y — > A"+^X. In local coordinates, 
we may write 

C{^) = L{x^,y^,v^)r+\ 

where L is a function on J^Y to which we also refer as the Lagrangian. 

The corresponding covariant Legendre transformation for a given Lagrangian den- 
sity JC : J^Y — > A"+^X is a fiber preserving map ¥C : J^Y — > J^Y* over Y, which is given 
by the first order vertical Taylor approximation to C: 



(F£(7),7') = >C(7) + ^ 



£(^ + e(7-7')), 



e=0 



where 7,7' e J^Y. In coordinates, F£(x^,y^, v^) = (x^, where 

III. THE HAMILTON-PONTRYAGIN PRINCIPLE FOR FIELD 
THEORIES 

In this section, we introduce the Hamilton-Pontryagin action principle for classical field 
theories. Let us first define the Pontryagin bundle 

M := J^Y Xy Z, 

as the fibered product over Y of the jet bundle J^Y and its dual Z (recall that we identified 
J^Y* with Z). This bundle plays a similar role as the Pontryagin bundle TQ © T*Q over a 
configuration manifold Q for the case of classical mechanics. Note that M is a bundle over 
Y and hence also over X. 

A. The Generahzed Energy Density 

Let jC : J^Y — > A"+^X be a Lagrangian density which is possibly degenerate. We define 
the generalized energy density associated to C by the map £ : M — >■ A'^'^^X defined as 



S (7, z) := ^(7, z)ri = {z, 7) - £(7) 
9 
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for (7,2;) e M. In local coordinates {x'^,y'^,v^,p,jf^ on M, the generalized energy density 
is represented by 



where E = p + p'^v^ — L{x^, y^, v^) is called the generalized energy on M. 

B. Pre-Multisymplectic Forms on M 

Recall the definition of the canonical forms © and Q := — d© on the dual jet bundle 
Z = J^y*, and let t^zm : M — > Z be the projection onto the second factor. The forms © 
and Q can be pulled back along t^zm to yield corresponding forms on M: 



Note that Q^m cannot be multisymplectic since it has a non-trivial kernel: for all v e TJ^Y 
we have that i^ilM = 0. As a result, we will refer to as the pre-multisymplectic 
{n + 2) -form on M. 

Finally, for any Lagrangian density C with associated energy density £, we introduce 
another pre-multisymplectic (n + 2)-form on M by 



C. The Hamilton-Pontryagin Principle 

Using the pre-multisymplectic forms Qm and Qg, we establish the Hamilton-Pontryagin 
variational principle for classical field theories. This principle is similar to the expression 
(1) for mechanical systems outlined in the introduction. 

Definition III.l. Consider a Lagrangian density C with associated generalized energy den- 
sity S. The Hamilton-Pontryagin action functional is defined as 



Jx 

where ip is a section of^xM M ^ X. In other words, ip can be written as ip — (7, z) where 
^ is a section o/7rx,jiy : J^Y — >■ X and z is a section of nxz Z ^ X. 



S^E{x^,y'',v^,p,p^^)d-+'x 
^ {p + P^i - L(x'^,y^v^)) d^+\ 



(14) 



©M T^zM^ and Qm '■— t^zm^- 




(15) 
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We have defined the Hamilton-Pontryagin principle in terms of the 1)- form ©m- By 
using the definition (10) and (13), we can rewrite the Hamilton-Pontryagin action functional 
as 



Jx 



where is a section of ttxy as before. An expression for S in local coordinates will be given 
below in Section HID. 

A vertical variation of a section -0 of ttxm : M — >■ X is a one-parameter family 
ipx '■ X ^ M oi sections, such that ipo — ip, and ipx agrees with ip outside a compact 
set U oi X. Such a variation can be constructed by considering a one-parameter family 
of diffeomorphisms ux : M ^ M such that uq is the identity, ux is the identity outside of 
\f{U) and ux preserves the fibration ttxm : M — )■ X, or in other words ttxm ° J^x = t^xm- 



'^XM 

Then, the composition ■i/'A = ^'a o is a vertical variation of ip. 

At the infinitesimal level, an infinitesimal variation of ^0 is a Tr^M-vertical vector field 
Vm : Im ^ — )■ VM, defined on the image of '4> in M. Here VM is the vertical subbundle of 
TM, defined by 

VM = {w e TM I Tt:xm{w) = 0} . 

We denote by X^(M) the module of the tt^m- vertical vector fields on M. 

We say that is a critical point of the action (15) if SS{ip) = for all variations of ip, 
where 



- _1 

~ dX 



X=0JU 

where U is a. compact subset of X that depends on the variation ipx- 

Proposition III. 2. A section ip of M is a critical point of the Hamilton-Pontryagin action 
functional if and only if ip satisfies the implicit Euler-Lagrange equations 



^*{VmJ ^s) = 0, for any Vm e X^(M). W 
Proof. Recall that a variation of ip is given hy tpx — i^x° where i/a : M — >■ M is the fiow 
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A. 
dX 



A=0 



A=0 JU 



of a vertical vector field Vm- Then, it follows that 

= [ r£vMi.^M-s) 

Ju 

= - / r{VM^{^M+d£))+ [ dir{VM-i{OM-s))) 

Ju Ju 

= - / r(VMj + f riVuJ (Om - s)) 

Ju JdU 

= - / riVM-i ns) 

Ju 

for all Vm, where we utilized Cartan's magic formula and Stokes' theorem and the fact that 
the TTjsfM-vertical vector field Vm : M — > VM is compactly supported inU C X. A standard 
argument then shows that t/j is a, critical point of 5" if and only if (16) holds. □ 

D. Coordinate Expressions 

Employing local coordinates {x^,y'^,v^,p,p'^) on M, the action functional is denoted by 

Ju 

^ X {^^ " <) + ^(-^ <) } d""^'-- 

A small coordinate calculation, similar to the one in the introduction, then shows that the 
implicit Euler- Lagrange equations (16) are locally given by 

dy^^^A d^^dL_ dL_ 

dxi^ ^M' dx,. dyA^ Pa q^a- y"-') 

E. Generalized Energy Constraint 

The coordinates on the Pontryagin bundle are {x^,y'^,v'l^,p,p'^), and the Hamilton- 
Pontryagin principle gives equations for the fields y^, and p^, but leaves p undetermined. 
In order to fix p, we impose the generalized energy constraint 

£; = p + p>^-L(a;^/,O = 0. 
12 
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Together with the third equation in (17), this naturally recovers the covariant Legendre 
transformation (12). We summarize the results up to this point in the following theorem. 

Theorem III. 3. The following statements for a section ip : X ^ M of ttxm satisfying 
E{ip) — are equivalent: 

(1) ip is a critical point of the Hamilton-Pontryagin action functional (15). 

(2) ■0*(Vm-I ^e) — for all irxM-vertical vector fields Vm on M. 

(3) ip satisfies the implicit Euler- Lagrange equations (17) together with the covariant Leg- 
endre transformation (12). 

F. Implicit Field Equations Using Multivector Fields 

We now rewrite the field equations (17) in a slightly different way, which will make the 
link with multi-Dirac structures in Section IV clearer. Up to this point, we have viewed 
the field equations as coming from the Hamilton-Pontryagin variational principle. We now 
show how these equations can be understood directly in terms of the multisymplectic form. 
A quick introduction to the geometry of multivector fields can be found in appendix A. 

To this end, we will use multivector fields on M (see Refs. 17, 18, and 30), which can 
be seen as the infinitesimal counterparts of sections i/j : X ^ M oi the Pontryagin bundle, 
in the sense that a multivector field gives rise to a set of first-order PDEs, whose integral 
sections are sections of the Pontryagin bundle. For the purpose of this paper, we will restrict 
our attention to multivector fields that are locally of the form 



component of the multivector field in the t'^-direction is zero. We will refer to multivector 
fields X which are locally given by (18) as partial multivector fields. 

The characterization (18) can be recast into an intrinsic form by noting that a partial 
multivector field A" is a map from M into f\^^^{TZ) with the following properties: 

• A" is decomposable: in other words, there exist n + 1 vector fields A"^ on M so that X 
can be written as the wedge product X — A^^]^A^. 




(18) 



where C^, C^^ and are local component functions of X on M. Note in particular that the 
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• A" is normalized, in the sense that ixiT^xMV) — 1 with rj the volume form on X. As a 
result, the coefficient of d/dx^ in (18) is unity. 

Since TZ is a subbundle of TM — T{J^Y) XtyTZ, we can identify a partial multivector 
field X with a regular multivector field on M which is J^Y -vertical, i.e. which has no 
component along d/dv^. More information about multivector fields and their properties 
can be found in Refs. 16 and 17. 

A section ip : X ^ M, given in local coordinates by 

i/j{x) = {x'',y^{x),y^{x),p{x),p'X{x)), 
is an integral section of X if it satisfies the following system of PDEs: 

which express the condition that the image of ip is tangent to X. We now claim that the 
implicit field equations (17) can be expressed using multivector fields and the generalized 
energy (14) as follows: 

Theorem III. 4. Let E he the generalized energy (14) and consider a multivector field X 
on M of the form (18) such that 

ix^M = (-i)"+M£;, (20) 

where dimX — n-\-l, and let ip : X ^ M be an integral section of X. In local coordinates, 
where ip{x) — {x^,y^{x),y^{x),p{x),p'^{x)), we then have that ip satisfies the following set 
of local equations: 

dp^^dL_ dy^^^ dL_ 

dxf^ dy^' dxi^ dvf ^ ' 

together with 

^,{p + PaV^-L)-^. (22) 

Proof. To show this, we compute the differential of the generahzed energy (14): 
, ^ 9^ , „ BE , , dE . BE , . dE , 

dL 



dxi^ 



+ (-^) + {p'a - dvi + v^dp^^ + dp. 
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and the contraction of the multivector field X given in (18) with the multisymplectic form 

= (-1)"+^ [{C^C'^x - C^CX - C,)dx^ + C^,dv\ + C\/y^ + dv\ . (23) 
As a result, we have that (20) holds if and only if 



dL 

Cax - C^^ -Cfj,^ (24) 



as well as 

= Ca,-§^, Pa-§^- (25) 

Now, let ip be an integral section of X, so that locally (19) holds. From (25), it follows that 
the component functions of ip satisfy 

and together with the third relation in (25) these equations are just the Euler-Lagrange 
equations (21) in implicit form. The expression (24) can now be rewritten by substituting 
into it the expression for from (25) and rewriting the second term using the relation 
— to obtain 



I^QyA A A^^ M g^^,- 

We now substitute the expressions from (19) determining the component functions of ip. 
The result is 



Using the Leibniz rule for the second term and regrouping the resulting terms, we then get 

dxi" dxi" dy^ dxi" dv^ dxi" ^^^^ ' ' 
which is precisely (22). □ 



IV. MULTI-DIRAC STRUCTURES AND FIELD THEORIES 

In this section, we show how the imphcit field equations '\x^m — (— l)""*"^d£' of theo- 
rem III. 4 can be described in terms of multi-Dirac structures, which are a graded analogue 
of the concept of Dirac structures introduced by Courant in Ref. 14. In order to keep the 
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technical side of the paper to a minimum, we focus on the applications of multi-Dirac struc- 
tures to field theories, and we leave proofs and a more thorough discussion to Ref. 38. We 
also point out that there are alternative approaches to the construction of a "field-theoretic 
analogue" of Dirac structures; see for instance Refs. 1 and 45. 

A. Pairings between Multivectors and Forms 

While multi-Dirac structures can in principle be defined on an arbitrary manifold M, we 
will focus for now on multi-Dirac structures on the Pontryagin bundle M ^ X. Recall that 
the dimension of the base space X is written as n -|- 1. 

Consider the spaces /\\TM) of i-multivector fields on M and l\^{T*M) of fc-forms on M. 
For k > I, there is a natural pairing between elements of S e f\\T*M) and X e /\\TM) 
given by 

(E, X) := ix^ e f\^-\T*M). (26) 
We now introduce the graded Pontryagin bundle of degree r over M as follows: 

where r = 1, n + 1. 

Using the pairing in equation (26), let us define the following antisymmetric and symmet- 
ric pairings between the elements of Pr and Pg as follows. For (A", E) e Pr and {X, E) e P,, 
where r, s = 1, . . . , n -|- 1, we put 

{{{X, E), {X, E)»_ ^ (i^E - (-l)"i^E) 

and 

{{{X, E), {X, E)»^ := i (i;pE + (-l)-i;,E) , 

each of which takes values in f\^^'^~'^''^{T*M). Hence, both of these pairing are identically 
zero whenever n -|- 2 < r -|- s. 

Let Vs he a, subbundle of Pg. The r-orthogonal complementary subbundle of Vg 
associated to the antisymmetric paring ((, ))_ is the subbundle {Vg)^'^ of Pr defined by 

(K)^'' = {(A',E)eP, |(((A',E),(A',S)»_ = for all {X,E)eVg}. 
Note that (Vg)-^'^ is a subbundle of Pr, and that (Fs)"*"'^ = Pr whenever n + 2 < r + s. 
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B. Almost Multi-Dirac Structures on Manifolds 

The definition of an almost multi-Dirac structure on M mimics the standard definition 
of Ref. 14 of Dirac structures. 

Definition IV. 1. An almost multi-Dirac structure of degree n-\-l on M is a direct 
sum Di® ■ ■ ■ ® Dn+i of subbundles 

Dr C Pr (for r = 1, . . . , n + 1^, 

which is (n + 1)-Lagrangian; namely, the sequence of the bundles Dj. are maximally 
{n + 1) -isotropic: 

Dr = (i^.)^'^ (27) 
for r,s — 1, n + 1 and where r + s < n + 2. 

When no confusion can arise, we will refer to the direct sum £)i © • • • © -D„+i simply as D. 
In the companion paper Ref. 38, we define a graded version of the Courant bracket, and we 
refer to integrable multi-Dirac structures as almost multi-Dirac structures which are closed 
under the graded Courant bracket. 

For the case of classical field theories, it will turn out that only the orthogonal comple- 
ments of the multi-Dirac structure of the "lowest" and "highest" order (namely, r — 1 and 
r — n-\- 1 respectively) play an essential role to formulate the field equations. 

Proposition IV. 2. Let D — Di ® ■ ■ ■ ® Dn+i be a multi-Dirac structure of degree n-\- 1. 
For any {X, S) G Dr and {X, S) G Dg, the following relation holds: 

i;eS - (-l)"i;,E = 0. (28) 
Proof. This is clear from the r-Lagrangian (maximally r-isotropic) property oi Dr- □ 

The above r-Lagrangian property of the multi-Dirac structure Dr in equation (28) may 
be understood as the field-theoretic analogue of the virtual work principle in mechanics 
and is related to Tellegen's theorem in electric circuits. 
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Example. Let M be a manifold with a multi-Dirac structure Di of degree 1. The isotropy 
property then becomes 

where Pi = TM © T*M and the pairing is given by 

(((A',E),(;e,E)»_ = ^(i;eS + i;.£), 

for all {X, S), {X, S) G Pi. This is nothing but the definition of an (almost) Dirac structure 
developed in Ref. 14. Our concept of multi-Dirac structures also includes, and is in fact 
equivalent to, the so-called higher-order Dirac structures of Ref. 45. 

C. Multi-Dirac Structures Induced by Differential Forms 

Consider an (n + 2)-form Qm on a manifold M. We will show that the graph of (in 
the sense defined below) defines a multi-Dirac structure of degree n + 1. This example of 
a multi-Dirac structure will be used in our subsequent treatment of classical field theories, 
where Qm will be the canonical multisymplectic form, but for now Qm can be an arbitrary 
form. 

Proposition IV. 3. Let Qm be an arbitrary {n + 2)-form on M and define the following 
subbundles D^, . . . , Dn+i, where C Pr is given by 

Dr^{{X,T.)ePr\ ix^M^T.} (29) 

for r = l,...,n-|-l. Then D — D^® ■ • ■ ® D^+i is a multi-Dirac structure of degree n-\-l 
on M. 

Proof. To prove that D — Di® ■ ■ ■ ® D^+i is a multi-Dirac structure of degree n-\-l, we 
need to check the isotropy property (27), namely, Dj. — D^'^ for all r, s = 1, . . . , n -|- 1, with 
r -|- s < n -|- 1. 

Let us first show that C D^''' . Let {X , E) e and (A', S) e D,. Then, it follows 
that 

(((A',E),(^,E)»_ 

= i{i;eE + (-l)-+M;,E} 
= 0, 
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since ip^ix^M = (—^Y^^xix^M- Thus, C D^''^. 

Next, let us show that D^'"" C Dr. Let {X, S) e D^''' . By definition of Dj'\ 

i^E+(-ir+4;,E = 

for all {X, E) e D^, i.e. A' e A''(7'^) such that i;tQM = S. Then, it follows 

i;eE + (-l)"+4^S = i;ei^nM + (-l)"+4^E 

= i;,{(-iri^QM + (-ir+'S} 
= 

for all X e KiTZ) where r + s < n + 2. Therefore, one has 

Thus, L*^''' C Dr and it follows that D,, = L*^''', so that D = Di ® ■ ■ ■ ® Dn+i is a 
multi-Dirac structure of degree n + 1 on M. □ 

Local Characterization of Multi-Dirac Structures. In Ref. 45, Zambon has given a com- 
plete local characterization of multi-Dirac structures. The following theorem is a combina- 
tion of corollary 3.9 and lemma 4.3 in that paper. 

Theorem IV. 4. Let T be an arbitrary vector space and consider a form oo G A'^^^ ^* ^'^^ 
a subspace S G T such that dim 5* < dimT — n or S = T. Then the sets 

Dr = {{X,E)\ X eSA (A'-^T), E - ixo; e A"+^-''5°}, 

r — 1, . . . ,n, form a multi-Dirac structure. Conversely, each multi-Dirac structure on T is 
of this form, for a suitable choice of S and ui. 

D. Applications to Implicit Field Theories 

We now consider Lagrangian field theories in the context of multi-Dirac structures. Con- 
sider a fiber bundle TTxr : Y ^ X and let M = J^Y Xy Z he the Pontryagin bundle. 
Let £ be a Lagrangian density on J^Y and denote the corresponding generalized energy 
density by £. Let ft be the canonical multisymplectic structure on Z and let be the 
pre-multisymplectic (n -|- 2)-form defined on the Pontryagin bundle by Qm = t^zm^- 
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Using Qjvf we can define a multi-Dirac structure D = Di © • • • ® D^+i of degree n + 1 on 
M by the construction of proposition IV.3. We refer to D as the canonical multi-Dirac 
structure on M. Exphcitly, D — Di ® • ■ ■ ® -Dn+i is given by 

Dr ^ {{X,i;,nM) \ X e /\'{TM)}, 

for r = 1, . . . , n + 1. As dD,M — 0, this multi-Dirac structure turns out be integrable; see 
Ref. 38. 

Prom theorem III. 4, we have that the implicit field equations can be written in terms of 
a partial multivector field X and the generalized energy E as 

This can be expressed in terms of the multi-Dirac structure induced by JIm by the condition 
{X, (— l)'^+^d£') G Dn+i- Note that the other components Di, . . . , Dn of the multi-Dirac 
structure are not needed to write down the field equations. 

Taking our cue from this observation, we now define a special class of field theories, 
whose field equations are specified in terms of a Lagrangian C, a multivector field X, and a 
multi-Dirac structure D on the Pontryagin bundle M. 

Definition IV. 5. Let D — Di ® ■ ■ ■ ® D^+i be a multi-Dirac structure of degree n + 1 on 
the Pontryagin bundle M and consider a Lagrangian density L with associated generalized 
energy E given by (14). A Lagrange- Dirac system for field theories is defined by a 
triple {X,E,Dn+i), where X is a normalized, decomposable partial vector field, so that 

(x,{-ir+'dE)eDn+i. 

Note that in the definition of a Lagrange- Dirac system, only the highest-order component 
Dn+i of the multi-Dirac structure appears. 

Theorem IV.6. For any Lagrange- Dirac system {X, E, Dn+i); the condition of conserva- 
tion of the generalized energy 

XldE^O (30) 

is satisfied. In other words, the generalized energy E is constant along the solutions of the 
Lagrange- Dirac system {X , E, Dn+i) (i.e. the integral manifolds of X). 
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Proof. The proof rehes on lemma A. 2 for decomposable vector fields. Using the field 
equations, the inner product on the left-hand side of (30) can be written as 

A'Ld£; = A'L(A'JQm) =0, 

according to Corollary A. 3. □ 

To see why the previous theorem implies energy conservation, decompose the multivector 
field X as in (18). Using lemma A. 2, the interior product can then be written as 

k 

A'Ld£; = ^(-1)^^+1 {X^,dE)X^ 
where A"^ is the n-multivector field obtained by deleting X^ from X, i.e., 

n+l 

= A 

A=l 

Since the multivector fields X^ are linearly independent, the energy conservation equation 
(30) then implies that {X^, dE) = for = 1, . . . , n + 1. In other words, the function E is 
constant along the integral manifolds of X. 

V. EXAMPLES 

In this section, we shall present some examples of Lagrange- Dirac field theories in the 
context of multi-Dirac structures. Namely, we show that the two examples of the scalar 
nonlinear wave equation and electromagnetism can be described as standard Lagrange-Dirac 
field theories. Furthermore, as an example of the Lagrange-Dirac field theories with nonholo- 
nomic constraints, we consider time- dependent mechanical systems with affine constraints. 
Lastly, we show that the Hu-Washizu variational principle of linear elastostatics can be 
understood in the context of the Hamilton-Pontryagin principle. 

A. Nonlinear Wave Equations 

Let us consider the scalar nonlinear wave equation in two dimensions as discussed in 
Ref. 5. The configuration bundle is t^xy :Y ^ X, with X = R^, so that n — dimX — 1 = 1, 
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and y = X X R. Local coordinates are given by {t,x) for [/ (Z X and {t,x,(f)) for Y, 
and hence {t,x,(t),Vt,Vx) for J^Y and {t,x,(f),p,p^,p^) for Z = J^Y*. The nonlinear wave 
equation is given by 

where is a section of ttxy and y : y — >■ M is a nonlinear potential. The Lagrangian for 
this equation is given by 

L{t, X, (p, vt, v^c) = ^ {vt - vl) + V{(j)) 

and the Lagrangian density is of course given by C — Ldt A dx. The generalized energy is 
then E + p% + p^^v^ - L{t, x, 0, Vt, v^). 

1. Hamilton-Pontryagin Principle 

The Hamilton-Pontryagin action functional (15) for the scalar field is given by 



Q K - + y{<t>) +P\<t>,t - Vt) - dt A dx. 



Taking variations of (f),Vt,Vx,p^,p^ (with the condition that d4> vanishes on the boundary of 
U), we obtain the implicit field equations 

d(f) <90 i ^. dp^ dp"" 

and imposing the generalized energy constraints E — 0, one has 

p = L{t,x,(f),vt,v^) -p*vt-p''v^. (32) 
The equations (31) are precisely the covariant equations obtained in Ref. 5. 

2. Multi-Dirac Structures 

The canonical pre-multisymplectic 3-form on M is given by 

flM = -c?P* Adcf) Adx + dp^ A d(f) A dt - dp A dt A dx, 

and one can define a multi-Dirac structure D = DiQ) D2 of degree 2 associated with Qm on 
M by the standard construction, where we set n = 1. We omit the definition of Di, since 
only D2 is important for the formulation of the dynamics: 

D2 = {{X, E) e P2 = A\TM) © T*M I ix^M = S} . 

22 



The Hamilton-Pontryagin Principle and Multi-Dirac Structures 



Here, A" is a partial multivector field locally given by 



d ^ d d d 

+ c„— + c:,— + c, 



where the variables x^, // = 0, 1, are given hj x'^ — t and x^ — x. The field equations can 
now be expressed as 

{X, -dE) e ^2, 

or in other words ix^M — —dE. In coordinates, these equations become 

- r „ - r - rf> - - 

acp ovq ovi 

(9L° dL 
Co + CqCi — CqCq = , ci + CqCi — CqCi = 'q^-> 

where 

r^v ^ dp_ 

^ dxi^' dxf^' dxf^' 

Hence, by simple rearrangements, we obtain the implicit Euler-Lagrange equation (31) for 
the nonlinear wave equation, together with 

^P ^ I T { 1 1 \ I \ 

By imposing the generahzed energy constraint = 0, we again obtain (32). 



B. Electromagnetism 

The multisymplectic description of electromagnetism can be found, among others, in 
Ref. 22. Here, we highlight the role of the Hamilton-Pontryagin variational principle and the 
associated multi-Dirac structure. The electromagnetic potential A — A^dx^ is a section of 
the bundle Y — T*X of one-forms on spacetime X. For the sake of simplicity, we let X be 
with the Minkowski metric, but curved spacetimes can be treated equally well. The bundle 
Y has coordinates {x^^,A^) while J^Y has coordinates {x^,A^,Ai^^^). The electromagnetic 
Lagrangian density is given by 

C{A,j'A) = ~dAA*dA = 

where the Hodge star operator * is associated to the Minkowski metric and where F^i, is the 
anti-symmetrization F^^, — ^4^ — Aj^^^. 
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1. Hamilton-Pontryagin Principle 

The Hamilton-Pontryagin action principle is given in coordinates by 



S 



u 



( ^ _ ^^^^^ _ If.„^^f,^ 



d X, 



where U is an open subset of X. The Hamilton-Pontryagin equations are given by 

and by ehminating p'^'*' lead to Maxwell's equations: d^F^" — 0. 

2. Multi-Dirac Structures 

The canonical multisymplectic form for electromagnetism is given in coordinates by 

= dA,, A dp^'" A d^'x^ -dp A d^'+^x. 

Let D = Di ® ■ ■ ■ ® D4 he the multi-Dirac structure induced by this form and consider in 
particular the component D4 of the highest degree. In terms of this bundle, the implicit 
Euler-Lagrange equations can be described by 

{X, -dE) e D4, 

where A" is a partial vector field of degree 4 given locally by 

4 



/ \ \QxI^ 

/i=i ^ 



d d d d 



dxi' ""^dA^ " dp''''' '"dp 
and E is the generalized energy density given by 

A„ = p+ p^^'A,,, + \f^.f^\ 

As before, a coordinate computation then shows that the coefficients of X are given by 

8T BT 
C ^ A C^" = — = v^'" = — = F^" 



and 
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where 



V 




and C., — 



dp 



UfJ, 



A similar computation as in the proof of Theorem III. 4 shows that these equations are 
equivalent to Maxwell's equations (33) in implicit form, together with the energy constraint 



C. Time-Dependent Mechanical Systems with Afflne Constraints 

We now show how time-dependent mechanics with linear nonholonomic constraints can 
be interpreted as a special case of Lagrange-Dirac field theories. Let Q be the configuration 
space of a mechanical system and consider the bundle tt^y : Y ^ X, where X = R 
(coordinatized by the time t) and y = R x Q, with coordinates (t, g*). The first jet bundle 
is given by J^Y = M x TQ, with local coordinates {t,q\v'). Let C : J^Y ^^{X) be a 
Lagrangian density and put C{j) = L{t,q\v^)dt. The dual jet bundle Z is a bundle over 
Y with local coordinates {t,q\pt,Pi) and can be identified with T*Y ^ T*R x T*Q. The 
canonical one-form on Z is given by © = Pidq^ + ptdt and the canonical two-form is given 
by Q = -d^ = dg' A dpi - dpt A dt. 

1. Non-Holonomic Affine Constraints 

Consider now a distribution A on y which is weakly horizontal, i.e., there exists a vertical 
bundle W (ZVY such that A®W ^TY (see Ref. 25). In local coordinates, let 



with a — 1, . . . ,k, be a basis for the annihilator A°. The distribution A gives rise to an 
affine subbundle of J^Y, denoted by C and given in coordinates by 



Intrinsically, C can be defined as follows. Recall that a jet 7 e J^Y can be viewed as a map 
from TjR to T(^t,q)y which is locally of the form ^ ^ dt ® {d/dt + v'd/dq') (see (9)). With 
this interpretation, C is the set of all 7 e J^Y such that Im7 C A. 



(22). 



:^Af{t,q)dq' + B'^{t,q)dt 



(34) 



C = {{t, q\ v') e J^Y I A'^it, qy + B^{t, q) = 0}. 
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The distribution A can now be hfted to the Pontryagin bundle M — J^Y © Z as follows. 
First, let Try jiy : J^Y — > y be the canonical jet bundle projection, given in coordinates by 

TTy^jiy (t, g*, w*) = {t,q''), and define now 

This is a distribution defined on the whole of J^Y, but only its values along C will be needed. 
We therefore denote the restriction of this bundle to C by A^. If the annihilator of A is given 
by (34), then an easy coordinate computation shows that A^ is spanned by the k contact 
forms 

r:^At{t,q){dq'-v'dt), (35) 

with a = 1, . . . ,k. Finally, we let A^ be the distribution defined along points of C x J^Y* C 
M, given by Am ^''''jiVm(^c), where njiy^M '■ M J^Y. The annihilator A^ is spanned 
by the puUback to M of the one- forms in (35). This is but one way of constructing A^; 
other approaches (which, unlike our construction, can be used for nonlinear constraints as 
well) can be found in Refs. 7 and 28, and the references therein. 



2. Dirac Structures 

We now construct a Dirac structure on the Pontryagin bundle M as follows: 

DM = {{X,^)eTM(BT*M\ i:-ixnMeAl, for X e Am} , 

where i^M '■— '^zm^- can be shown directly that Dm is a Dirac structure (see Ref. 41); 
alternatively, note that Dm is of the form described in the theorem IV.4 and hence auto- 
matically determines a (multi-)Dirac structure. 

Now, let A* : M ^ TZ be a partial vector field as in (18). In local coordinates, we have 

d . d . d . d 

^ + 5 ^ + Pi-^ + Pt-^- 
at oq^ opi opt 

On the other hand, the generahzed energy E is given by £^ = Edt, where 

E{t, q\ v\pt,Pi) =pt+ PiV' - Lit, q\ v') 
and the differential of E is locally given by 

- — j dt + f J dq' + I Pi - — ] dv' + dpt + v'dpi. 
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Thus, the Lagrange-Dirac equations for nonholonomic time-dependent mechanics are 
given by the triple (A", E, D) which satisfies the condition 

{X, dE) e D. (36) 

This is in turn equivalent to the intrinsic nonholonomic Lagrange-Dirac equations: 

ix^ - d£; e and X e Am- 

In coordinates, one obtains 

^-v\ p,^- + KA:, p^ = —, (37) 
together with the affine constraints Al{t, q)v^ + B^{t, q) — Q and the energy conservation law 

By imposing the constraint that = 0, we find that pt = L — piv\ Note that the equations 
of motion (37) are the analogues of the Dirac equations (4) for affine constraints. 

One can apply a similar approach to construct multi- Dirac structures for classical field 
theories with nonholonomic constraints. After constructing the distribution A^ in a suit- 
able way (see e.g. Ref. 37), the nonholonomic multi-Dirac structure will then be given 
again by theorem IV.4 and the nonholonomic field equations will be given by (36). At this 
point, we emphasize that in the current description of nonholonomic multi-Dirac structures 
the constraints are assumed to be linear/ affine in the multi- velocities, while A^ can be 
constructed for nonlinear constraints as well. Most physically relevant nonholonomic field 
theories, however, have nonlinear constraints, and the study of nonlinear nonholonomic field 
theories using multi-Dirac structures will be the subject of future work. 



D. Elastostatics 

To conclude, we show how the Hu-Washizu principle from linear elastostatics (see 
Ref. 40) can be understood as a Hamilton-Pontryagin variational principle. Let B be an open 
subset of representing an elastic medium and consider the bundle tt : Y :— B x M.^ ^ B. 
We denote coordinates on J^Y by {X^ , (f)\ F}), where the fiber coordinates F} represent the 
deformation gradient of the material. The dual bundle Z has coordinates {X^, (t>^,P, t^I), 
where tt/ will be seen to represent the first Piola-Kirchhoff tensor. 
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Furthermore, we let : J^y — > R be the stored- energy function of the material, 
and we assume that 1^ is a function of F} only (this assumption can easily be relaxed) . We 
denote the external body forces acting on the medium by bi{X) and we assume that there 

are traction boundary conditions acting on the boundary dB, given by n-nj = Ti, where 
nj is the unit normal to the boundary and Tj(X) is a prescribed traction vector. Dirichlet 
boundary conditions can easily be incorporated by restricting the space of fields as in Ref. 31, 
though we will not do so here. 



1. Hamilton-Pontryagin Principle 



Taking W as our Lagrangian, the Hamilton-Pontryagin principle becomes 



S{<l>,F,n) 



F})+W{F)-k{X)<j>\X) 



dV 



T,{X)(t>\X)dS, (38) 



dB 



where we have incorporated the traction boundary conditions through the surface integral. 
Due to the presence of the boundary integral, this action functional is more general than 
the Hamilton-Pontryagin principle that we dealt with in Section HI, but for the purpose 
of taking variations it can be treated in exactly the same way. In Section V D 3, we will 
formulate a different way of incorporating the boundary conditions, which can easily be 
generalized to the case of elastostatics on a nonlinear manifold. 
Taking variations of (38) with respect to 0, F, and tt results in 



5S^ 



dX' 



Fl 54 + 



dW 



+ 



OF} 

I {T:lni-Ti{X))8(t>\QsdS, 

JdB 



SF} 



dn: 



dX 



dV 



which leads to the following equations: 



dX' 



dXi' 



and tt/ 



dW 



together with the boundary conditions 



(39) 



(40) 



TTj Hi — Ti (on dB) , 

where ni is the outward normal vector to dB. The last equation in (39) shows that tt/ 
can be identified with the first Piola-Kirchhoff tensor, the second equation is nothing but 
the definition of the deformation gradient, while the first equation is the Navier-Cauchy 
equation. 
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2. Hu-Washizu Principle 

The variational principle (38) is not yet the Hu-Washizu variational principle. To obtain 
the latter, we work in the context of linear elasticity. Write the elastic field 0' as (f)^{x) — 
x'^ + u^{x), where u = (m^, w^, u^) is the displacement field, and where material and spatial 
points are both denoted by x\ For the stored energy function, we take a quadratic expression 
of the form 

W{e)^lciie}et with e} = ^ (^^^ + x.^) , 

where are referred to as the strain tensor, while the constitutive tensor cH character- 
izes the elastic medium. We now substitute these expressions into (38). Up to a global sign 
and ignoring boundary conditions and external forces for the sake of simplicity, we obtain 
the Hu-Washizu principle for linear elastostatics: 



S{u,e,p) = / 
Jb 



dV. 



Here, we have used the fact that the stress tensor pj = dW{e)/de^j is symmetric in both 
indices. 

It goes without saying that this derivation relies heavily on the Euclidian structure of the 
space in which the body moves. Nevertheless, it is possible to derive a Hu-Washizu principle 
for nonlinear elasticity in terms of the Green tensor which is valid on arbitrary Riemannian 
manifolds. This modified principle can be obtained from the Hamilton-Pontryagin principle 
by means of symmetry reduction and will be the subject of a forthcoming paper. 



3. The Hamilton-Pontryagin Principle with Boundary Conditions 

To finish, we wish to point out that the boundary term in (38) only makes sense when 
the space in which the medium moves is Euclidian. When this is not the case, the boundary 

conditions (40) can be incorporated in a different, more intrinsic way which is reminiscent 
of the Lagrangc-d'Alembcrt principle in mechanics (see for instance Ref. 26). We first let S 
be the action integral (38), but without the boundary term: 

5(0, F, tt) = / 
Jb 
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where for the sake of convenience we have also set the external body forces bi{X) equal to 
zero. We require now that the physical fields (0, F, tt) are those for which S satisfies 



for arbitrary variations {6(f), 6F, Sn). This variational principle reproduces the field equations 
(39) and the boundary conditions (40), but has the advantage of being closer in spirit to 
the Hamilton-Pontryagin principle, and of being more intrinsic: when the elastic fields take 
value in a manifold Q, the boundary term in (38) is only locally defined. By contrast, (42) 
holds globally true even when the target space is a manifold: the variations 50 are then 
elements of TQ, while the traction vector r can be seen as an element of T*Q, and the 
pairing between both in (42) is the standard duality pairing. 

We now show how the field equations (39), together with the boundary conditions (40), 
can be expressed intrinsically in terms of the multisymplectic form as in Theorem 111.3. To 
this end, let /3 be a map from the boundary dB to the dual jet bundle Z, given in coordinates 



which will encode the boundary conditions, in a form to be determined later. We recall that 
in this example the Pontryagin bundle M has coordinates {X^ , 4>\ F},p,tt^), and we denote 
a generic section of M by ^0. 

With S given by (41), consider now the modified Hamilton-Pontryagin principle 



for any finite variation ipx of ijj, with infinitesimal variation given by the vertical vector 
field Vm (see Section lllC, where this notation was introduced). Following the proof of 
Proposition III. 2, but keeping the boundary terms, we see that this variational principle is 
equivalent to 




(42) 



by 





Jb 

leading to the field equations 




and the boundary conditions 



'M 



(43) 
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for all vertical vector fields Vm- To see that these expressions reduce to (39) and (40), 
observe that ©m and E are given in the case of elastostatics by 

Om = 7r/#^ A 6^X1 + and E = (tt^ F}+p- W)(P+^X. 

It then follows easily that the field equations are given by (39). For the boundary conditions, 
note that the vertical vector fields are spanned by {d/d(f)\ d/dF}, d/dn- , d/dp}. When 
evaluated on these basis elements, (43) vanishes unless Vm = d/d(p\ In that case, we have 

= ri'^a/aA^M -E-P)) = r{{< - BDdTXj). 

If B is equipped with a Riemannian metric gu, the last expression can be rewritten by 
noting that ip*{d^Xi) — g~^/^nidS, where nj is the unit normal vector, g is the determinant 
of the metric, and dS is the induced volume form on dB. The boundary conditions (43) 
then become nfni — B(ni, which coincides with (40) if we choose /3 so that Bfni — r^. 

VI. CONCLUSIONS AND FUTURE WORK 

In this paper, we have developed the Hamilton-Pontryagin principle for field theories 
and we have studied the relation between the field equations obtained from the Hamilton- 
Pontryagin principle and Lagrange-Dirac field theories, formulated in terms of multi-Dirac 
structures. These ideas are shown to be a natural extension of Dirac structures and 
Lagrangc-Dirac dynamical systems in mechanics. Finally, we have given a few illustrative 
examples of Lagrange-Dirac field theories: the Klein-Gordon equation. Maxwell's equations, 
the case of time-dependent mechanics with affine nonholonomic constraints, and elastostat- 
ics. 

Multi-Dirac Structures for Field Theories with Boundary Conditions. We have shown in 
this paper how the Hamilton-Pontryagin principle can be augmented in order to deal with 
field theories with non-trivial (Neumann) boundary conditions. We are currently working 
on an extended concept of multi-Dirac structure, in which the boundary conditions are 
automatically incorporated. As we discuss below, we expect that our approach will be 
related to the Stokes-Dirac structures of van der Schaft and Maschke (Ref. 34) used in 
boundary control theory. 
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Space-Time Slicings and the Adjoint Formalism. The approach in this paper are covari- 
ant in the sense that no distinction is made between time and the spatial variables on the 
base space. In fact, as we have pointed out before, the configuration bundle ttxy :Y ^ X 
can be arbitrary and in particular X does not have to represent spacetime. While this 
approach has various advantages (see Ref. 5, 22, and 24 for an overview), not in the least 
that the resulting structures are all finite-dimensional, it is often necessary to single out one 
coordinate direction as time, for instance when considering the initial- value problem for field 
theories. 

In related work (see Ref. 11) we will show that when such a space-time slicing is chosen (as 
in Ref. 2 and 23), the Hamilton-Pontryagin variational principle yields the field equations in 
adjoint form., bringing out the distinction between the dynamical variables of the theory (the 
physical fields) and the kinematic variables, which are present due to gauge freedom. For 
instance, in the case of relativity, the former are the metric and the extrinsic curvature on 
a spatial hypersurface, the latter are the lapse and shift, and the adjoint formalism reduces 
to the familiar ADM equations of relativity. 

Furthermore, the choice of a space-time splitting induces a correspondence between multi- 
Dirac structures on the jet bundle and certain infinite-dimensional Dirac structures on the 
space of all fields. The latter are in turn related to the concept of Stokes-Dirac structures 
(Ref. 34). The precise link between Stokes-Dirac structures and infinite-dimensional Dirac 
structures involves Poisson reduction and was established in Ref. 39. 

Discretizations of Multi-Dirac Structures and Field Theories. Leok and Ohsawa (Ref. 27) 
recently proposed a concept of discrete Dirac structure in mechanics. At the same time, it 
has been shown in Ref. 4 that discrete versions of the Hamilton-Pontryagin principle can 
be used to derive very accurate geometric integrators for mechanical systems. We propose 
to explore the field-theoretic version of these discretizations. Since multi- Dirac structures 
and variational principles are a natural way to describe field theories even for degenerate 
Lagrangians, this approach is especially promising. In particular, we will focus on field 
theories whose covariance group includes the diffeomorphism group of the base space X, 
focusing on the link with the adjoint formalism. 

Symmetry and Reduction of Multi-Dirac Structures. The presence of non-trivial symme- 
try groups is an essential feature of many classical field theories, such as electromagnetism, 
Yang-Mills theories, and elasticity. More complex examples include the theory of perfect fiu- 
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ids and Einstein's theory of relativity. It is often useful to obtain a form of the field equations 
from which the symmetry has been factored out; this is for instance what happens in writing 
Maxwell's equations in terms of the field strength F^y rather than the vector potential A^. 
From a bundle-theoretic point of view, various approaches of symmetry reduction for field 
theories have been proposed (see Refs. 9, 10, and 19 and the references therein). We pro- 
pose to develop a reduction theory for multi-Dirac structures along the lines of the theory of 
symmetry reduction for Dirac structures in classical mechanics developed by Refs. 3, 43, and 
44. 
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Appendix A: Generalities about Multi- Vector Fields 
1. Definitions 

Let M be a manifold of dimension n. A /c-multivector field (where A; < n) on M is a 
section X of the A;- fold exterior power /\^{TM) of the tangent bundle. The module of all 
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/c-multivector fields is denoted by X'^{M). 

We define the interior product of multivector fields and forms as follows. Let X e 
X'^(M) and a e n^{M), where k < I. The left interior product X} a is then the unique 

(/ — /c)-form such that 

{XAa){X') = a{X ^X') 

for all {I — A;)-multivectors X' . Similarly, the mght interior product of a multivector 
X e X''{M) and a form ^ e Q"*(M), with k > m, is the unique {k — m)-multivector Xif3 
with the property that 

(A'L/3)J7 = (/5A7)(A') 

for all (k — m)-forms 7 on M. 

Further information about multivectors and interior products can be found in (for in- 
stance) Ref. 30 and 36. Prom these papers, we quote the following two results on the 
relation between the interior product and the wedge product of vector fields and forms, 
respectively. 

Lemma A.l (Ref. 36). Let X e X''{M) and a e Q^M). Then 

{X AX')la = {X I a) A X' + {-!)'' X A {X' I a) 
for any multivector X' , and 

X\ {a A 13) = (A'La)J/9 + (-l)'^a A (A'J/3) 

for any form (3. 

When the multi-vector field is decomposable, the right interior product with a one-form 
can easily be computed by means of the following formula, the proof of which follows from 
the previous theorem. 

Lemma A. 2. Let X be a k -multi-vector field which is decomposable in the sense that 
X = A^=i '^/i? where X^ (fi = 1, . . . , k) are vector fields and consider a one-form a. Then 

k 

where Xf^ is the {k — 1) -vector field obtained by deleting X^, i.e. 

n+l 
A=l 
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Corollary A. 3. Let X be a decomposable k-multivector and a an l-form, where I > k. Then 
Xl(X}a) = 0. 

Proof. Let (5 — Xja. With the notations of the previous lemma, we have that 

k 

but (A"^, /3) = for all = 1, . . . , k, since {X^, /3) — ix^i^n+i ^^a = 0. □ 

In this paper, we also let X^y (where [i ^ v) \>% the {k — 2)-multivector obtained by 
deleting both X^ and X^,: 

n+l 

A=l 

li II < u, we then have that 

Xx A x^. = i-ir+'s^xx, - i-iy+XxX^. 

U H > u, note that X^^ = X^f^. The above formula can then be applied. 
2. Proof of (23) 

The computation of equation (23) is somewhat involved and depends on a number of 
coordinate identities, listed here. First, recall that X is an (n + l)-multivector field with 
local expression (18), and that the multi-symplectic form Q is locally given by equation (11). 

Recall the volume form 77 on X is locally expressed by d'^'^^x and denote 

d d 

The following contractions will be useful: 

A'^Jd"x,= (-l)^+'^5^, (Al) 

and, for < i/, 

X^u^ d^xx = {-ir^''+''{Syxdx>' + S^xdx"). (A2) 

Now return to the calculation of equation (23) and the contraction of X with is now 
given by 

X} Qm = X-l (dy^ A dp'X A d"x^) - X} (dp A d"+^x), 
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and the two terms on the right-hand side will be calculated separately. 
We begin with the first term. Using lemma A.l, we have 

Both terms can be calculated using lemma A. 2 and using (Al) and (A2). After some 
rearrangements, the result is that 

A'J (dy^ A dp^^ A d%) = (-1)"+^ [{C^C\^ - CtC\^)dx^ + C^dp^^] . 

For the second term, we use the same techniques to conclude that 

A'J {dp A d''+^x) = (-1)"+^ [dp - Ci.dx''] . 

Putting the results for both terms together, we obtain (23). 
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